BRANCHING ANALYSIS OF A COUNTABLE FAMILY 
OF GLOBAL SIMILARITY SOLUTIONS 
OF A FOURTH-ORDER THIN FILM EQUATION 



P. ALVAREZ-CAUDEVILLA AND VICTOR A. GALAKTIONOV 

Abstract. The main goal of the paper is to justify that source-type and other global- 
in-time similarity solutions of the Cauchy problem for the fourth-order thin film equation 

(0.1) ut = -V • (M"VAu) in R N x K+, where n > 0, N>1, 

can be obtained by a continuous deformation (a homotopy path) as n — > + by reducing 
to similarity solutions (given by eigenfunctions of a rescaled linear operator B) of the 
classic bi-harmonic equation 

(0.2) u t = -A 2 u in R N x M + , where B = -A 2 + \ y ■ V + f I. 

This approach leads to a countable family of various global similarity patterns of (|0.1| 
and describes their oscillatory sign-changing behaviour by using the known asymptotic 
properties of the fundamental solution of (|0.2j) . The branching from n = + for (|0.1[) 
requires Hermitian spectral theory for a pair {B,B*} of non-self adjoint operators and 
leads to a number of difficult mathematical problems. These include, as a key part, the 
problem of multiplicity of solutions, which is under particular scrutiny. 



1. INTRODUCTION: TFES, CONNECTIONS WITH CLASSIC PDE THEORY, LAYOUT 

1.1. Main models, their applications, and preliminaries. We study the global-in- 
time behaviour of solutions of the fourth-order quasilinear evolution equation of parabolic 
type, called the thin film equation (the TFE-4, in short) 

(1.1) u t = -V ■ (\u\ n VAu) in R N x R + , 

where V = grad x , A = V • V stands for the Laplace operator in Mr, and n > is a real 
parameter. Fourth- and sixth-order TFEs (the TFE-6) having a similar form, 

(1.2) ut = V-(M n VA 2 u), 

as well as more complicated doubly nonlinear degenerate parabolic models (see typical 
examples in [29J), have various applications in thin film, lubrication theory, and in several 
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other hydro dynamic-type problems. We refer e.g., to [131 [151 [251 [26] f° r most recent sur- 
veys and for extended lists of references concerning physical derivations of various models, 
key mathematical results, and further applications. Since the 1980s, such equations also 
play a quite special role in nonlinear PDE theory to be discussed in greater detail below. 

The TFE-4 (11.11) is written for solutions of changing sign, which can occur in the 
Cauchy problem (the CP) and also in some free-boundary problems (FBPs); see proper 
settings shortly. It is worth mentioning that non-negative solutions with compact support 
of various FBPs are mostly physically relevant, and that the pioneering mathematical 
approaches by Bernis and Friedman in 1990 [5] were developed mainly for such solutions. 

However, solutions of changing sign have been already under scrutiny for a few years (see 
[SKEMIE]), which in particular can have some biological motivations [20], to say nothing of 
general PDE theory. It turned out that these classes of the so-called "oscillatory solutions 
of changing sign" of (II. ip were rather difficult to tackle rigorously by standard and classic 
methods. Moreover, even their self-similar (i.e., ODE) representatives can lead to several 
surprises in trying to describe sign-changing features close to interfaces; see [H] for a 
collection of such hard properties. It turned out also that, for better understanding of 
such singular oscillatory properties of solutions of the CP for (11.11) . it is fruitful to consider 
the (homotopic) limit n — > + , owing to Hermitian spectral theory developed in [12] for 
a pair {B, B*} of linear rescaled operators for n = 0, i.e., for the bi-harmonic equation 

(1.3) u t = -A 2 u in 1^x1+, B = -A 2 + ±yV + f /, B* = -A 2 - \ y ■ V, 

which will always be key for our further analysis. 

In the present paper, using this continuity/homotopy deformation approach "n — > + ", 
we shall focus our analysis to the Cauchy problem for (II. lj) for exponents n > 0, which 
are assumed to be sufficiently small. Some key and necessary references will be presented 
later on. We study the large time behaviour of the solutions of (11.11) . To this end, we will 
use some natural connections with a similar analysis for more complicated models such 
as the limit unstable fourth-order thin film equation (the unstable TFE-4): 

(1.4) ut = -V • (|u| n VAu) - A(H p -\i) , 

with the unstable homogeneous second-order diffusion term, where p > 1 is a fixed expo- 
nent; see [13] for physical motivations, references, and other basics. Here, (11.41) represents 
a fourth-order nonlinear parabolic equation with the backward (unstable) diffusion term 
in the second-order operator. Blow-up and global self-similar solutions of (11.41) have been 
extensively studied in [131 EHJ f° r the unstable TFE-4 (ll.4p and in [T5l [16] for the unstable 
TFE-6, 

(1.5) u t = V ■ {\u\ n VA 2 u) - AQul^u) , 

where further references and other related higher-order TFEs can be found. 

From the application point of view, it is well known (see references to surveys above) 
that (11.11) and (11.41) arises in numerous areas. In particular, those equations model the 
dynamics of a thin film of viscous fluid, as the spreading of a liquid film along a surface, 
where u stands the height of the film (then clearly u > that naturally leads to a FBP 
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setting). In particular, when n = 3 we are dealing with a problem in the context of 
lubrication theory for thin viscous films that are driven by surface tension and when 
n = 1 with Hele-Shaw flows. It is also important to note that, in ( \1A\\ . the fourth-order 
term reflects surface tension effects and the second-order term can reflect gravity, van der 
Waals interactions, thermocapillary effects, or geometry of the solid substrate. 

Finally, in order to summarize let us mention again that higher-order semilinear and 
quasilinear parabolic equations occur in applications to thin film theory, nonlinear diffu- 
sion, lubrication theory, flame and wave propagation (the Kuramoto-Sivashinsky equa- 
tion and the extended Fisher-Kolmogorov equation), phase transition at critical Lifshitz 
points and bi-stable systems (see Peletier-Troy [3S] for further details, models, and re- 
sults). Moreover, in the special situation when n = we should notice that (II. 4p is the 
well known unstable Cahn-Hilliard equation (the CHE) 

(1.6) u t = -A 2 u- A^u^u) in R N x R + ; 

see main references and full details in [T7] . The unstable semilinear model (11.61) and similar 
stable ones are much better known and are connected with several new applications that 
have increased the interest of the study of their quasilinear TFE extensions (11.41) and 
consequently of (II. ip . Note that, without any doubts, the semilinear CHE (11.61) in the CP 
setting, admits oscillatory solutions of changing sign, though with no finite (i.e., infinite) 
interfaces. As our main goal, we plan to extend those properties of the CHE (11.61) to the 
TFEs for small n > 0, where oscillations begin to concentrate at finite interfaces. 

1.2. A digression to reaction-diffusion theory. Furthermore, in the CP setting for 
(II .6p . one can write (11.61) in the form 

(1.7) Au t = Au + u p , where A : = (-A) -1 

is a standard positive operator, so that (II. 7p is a pseudo-parabolic second-order equa- 
tion. Such equations had been widely studied since the 1970's, and nowadays are well- 
understood with both the existence and uniqueness results of local and global classical 
or blow-up solutions obtained, first blow-up results for such pseudo-parabolic PDEs were 
due to Levine in 1973, [33] • As was noted in [17], there are some similarities between ( 11.71) 
and the classical semilinear heat equation from combustion theory (the solid fuel model) 

(1.8) u t = Au + u p in R N x R + , p > 1 , 

with N > 1, especially for the blow-up singularity formation phenomena. Mathemati- 
cal literature devoted to the study of (jl.8p include a huge number of papers published 
since Fujita's classic papers in the 1960s, and this remarkable history has been already 
explained in a dozen of well-known monographs; we refer to [TJ HH HU [3U [38], HO, HD H2], 
where further extensions and references can be found. In particular, concerning blow-up 
patterns, a complete description of all possible types of blow-up have been achieved for 
some ranges of the parameters p and N, especially in the subcritical Sobolev range 



Note that, for p > psoboiev; such a classification of blow-up patterns is far away from being 
complete, with a number difficult open problems posed. 

Nevertheless, using a standard pseudo-parabolic form ( 11. 7ft of the C-H flows could rise 
a hope to apply a huge experience achieved earlier for classic reaction-diffusion models 
fll.8p . though this is not expected to be that straightforward. 

1.3. A digression to porous medium equation: homotopy to the heat equation. 

Returning to the TFEs, note that the unstable nonlinear operator in (11.41) gives us the 
notorious classic porous medium equation (the PME-2, but posed backwards in time), 

(1.9) u t = A(|n| n_1 M) (for convenience, p is replaced by n as in (II. ip ) , 

which derives its name from the role in the description of flows in porous media. Parabolic 
PDE models in filtration theory of liquids and gases in porous media were derived by 
Leibenzon in the 1920s and 1930s, as well as by Richard's (1931), and Muskat (1937). In 
fact, modern filtration theory goes back to the beginning of the twentieth century initiated 
in the works by N.Ye. Zhukovskii, who is better known for his fundamental research in 
aerodynamics, hydrodynamics, and ODE theory (on his pioneering non-oscillation test 
in 1892, see [18, p. 19]). His contribution to "theory of ground waters" is explained in 
P.Ya. Kochina's paper [3TJ. For an extended list of references on this subject and more 
filtration history, see [T9] . 

It is well understood that, for any n > 1, (jl.9p has a family of exact self-similar 
compactly supported source type solutions (the ZKB ones from 1950s), which describe 
the large time behaviour of compactly supported solutions with conservation of mass in 
the case on a non-zero mass, i.e., 

/ u(x, t) dx 7^ 0. 

On the other hand, (11.91) also admits a countable (at least) family of other similarity 
solutions; see [21] for key references and most recent results. 

The PME-2 (ll.9p can be interpreted as a nonlinear degenerate version of the classic 
heat equation for n = 0, 

u t = Au in R N x R + . 

Note that passing to the limit n — > + in (11.91) for nonnegative solutions used to be 
a difficult mathematical problem in the 1970s-80s, which exhibited typical (but clearly 
simpler than in the TFE case) features of a "homotopy" transformation of PDEs. This 
study in ID was initiated by Kalashnikov in 1978 |28J . Further detailed results in were 
obtained in [2]; see also [10]. More recent involved estimates were obtained in [361 [37] 
for the ID PME-2 (II .9p establishing the rate of convergence of solutions as n — > , 
such as 0(n) as n — > 0~ (i.e, from n < 0, the fast diffusion range, where solutions are 
smoother) in L X (R) [35], and 0(n 2 ) as n -> 0+ in L 2 (R x (0,T)) [37J. However, the most 
of such convergence results are obtained for nonnegative solutions of (11.91) . For solutions 
of changing sign, there are some open problems; see [21] for references and further details. 

Note that, as customary, any kind of detailed asymptotic analysis for higher-order 
equations is much more difficult than for second-order counterparts in view of the lack of 
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Maximum Principle, comparison, and order preserving semigroups and potential proper- 
ties of the operators involved. Thus, practically all the existing methods for the PME-2 
1Q|) are not applicable to the TFE-4 flU)) or (Ojl . 

Thus, in the twenty-first century, higher-order TFEs such as (II. ip and (11.21) . though 
looking like a natural counterpart /extension of the PME-2 (11.91) . corresponding mathe- 
matical TFE theory is more complicated with several problems, remaining open still. 

1.4. Main approaches, results, and layout. It is worth mentioning that, unlike the 
FBPs, studied in hundreds of papers since 1980s (see [25J and [14J for key references and 
alternative versions of uniqueness approaches), thin film theory for the Cauchy problem 
for (11.11) or (ll.4p led recently to a number of difficult open problems and is not still fully 
developed; see the above references as a guide to main difficulties and ideas. In fact, 
the concept of proper solutions is still rather obscure for the Cauchy problem, since any 
classic or standard notions of weak- mild- generalized-... solutions fail in the CP setting. 

In this work, we perform a more systematic than before analysis of the behaviour of the 
similarity solutions through a so-called homotopic approach (branching from n — 0) via 
branching theory, using the Lyapunov-Schmidt methods for obtaining relevant results and 
properties for the solutions of the self-similar equation associated with (II. ip and, hence, 
for the proper solutions of (II. ip . Overall, loosely speaking, this approach is characterized 
as follows: good proper (similarity or not) solutions of the Cauchy problem for the TFE 
(II. ip are those that can be continuously deformed (via a homotopic path) as n — > + to 
the corresponding solutions of the bi-harmonic equation (ll.3p . which will play a crucial 
role in the subsequent analysis. This homotopic-like approach is based upon the spectral 
properties known for the linear counterpart (11.31) of the TFE (II. ip . Moreover, owing to 
the oscillatory character of the solutions of the bi-harmonic equation (11.31) being a "limit 
case" of the TFE (11.11) . close to the interfaces, this homotopy study exhibits a typical 
difficulty concerning the desired structure of the transversal zeros of solutions, at least 
for small n > 0. Proving such a transversality zero property is a difficult open problem, 
though qualitatively, this was rather well understood, [15] . 

Indeed, we ascertain through an analytic "homotopy", which is understood as just 
the existence of a continuity as n — > + , over the CP performed in the last section of 
this paper that the solutions of (11.11) are homotopic to the solutions of the bi-harmonic 
equation (II. 3p in a weak sense. However, we must admit that this does not solve the 
problem of uniqueness of solutions of the CP (see details in Section [5]), since the final 
identification of the solutions obtained via analytic e-regularization and passing to the 
limit e — > + remains not fully understood. Overall, it seems that e-regularizations 
of solutions of the TFEs via families of uniformly parabolic analytic flows, which was 
a powerful and successful tool for second-order degenerate parabolic PDEs (such as the 
PME-2 (II. 9p ). for TFEs again leads to difficult boundary layer-type problems that remain 
open in a sufficient generality. 

Some parts of the study of the thin film equation (jl.4p can be performed in similar 
lines, though a full homotopy approach would include the passage to the limit p — > 1 + , 
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leading to the limit linear equation (not treated here) 

u t = — A 2 u — Au. 

Thus, the layout of the paper is as follows: 

(I) Study of a countable family of global self-similar solutions of (11.11) via their branching 
from eigenspaces at n = + , Sections EHU and 

(II) Some general aspects of the CP for (II. ip by another homotopy approach, Section [5j 

2. Problem setting and self-similar solutions 

2.1. The FBP and CP. For both the FBP and the CP, the solutions are assumed to 
satisfy standard free-boundary conditions: 

{u = 0, zero-height, 
V« = 0, zero contact angle, 

— n ■ X7(\u\ n Au) = 0, conservation of mass (zero-flux) 

at the singularity surface (interface) T [u], which is the lateral boundary of 

supp mcR w xR + , TV > 1 , 

where n stands for the unit outward normal to To^]- Note that, for sufficiently smooth 
interfaces, the condition on the flux can be read as 

lim n- V(\u\ n Au) = 0. 
dist(x,roM)lo 

It is key that, for the CP, the solutions are assumed to be "smoother" at the interface than 
those for the FBP, i.e., (12. ID are not sufficient to define their regularity. These maximal 
regularity issues for the CP, leading to oscillatory solutions, are under scrutiny in [TJ] . 
Next, denote by 

M(t) := J u(x,t) dx 

the mass of the solution, where integration is performed over smooth support (M. N is 
allowed for the CP only). Then, differentiating M(t) with respect to t and applying 
the Divergence Theorem (under natural regularity assumptions on solutions and free 
boundary), we have that 

m i n . v(i«rA«) . 

r n{t} 

The mass is conserved if J(t) = 0, which is assured by the flux condition in (12.11) . 

The problem is completed with bounded, smooth, integrable, compactly supported 
initial data 

(2.2) u(x, 0) = u Q {x) in T [u] n {t = 0}. 

In the CP for (11.11) in M. N x R +) one needs to pose bounded compactly supported initial 
data (12. 2p prescribed in M. N . Then, under the same zero flux condition at finite interfaces 
(to be established separately), the mass is preserved. 
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2.2. Global similarity solutions: towards a nonlinear eigenvalue problem. We 

now begin to specify the self-similar solutions of the equation (11.11) . which are admitted 
due to its natural scaling-invariant nature. In the case of the mass being conserved, we 
have global in time source-type solutions. 
Using the following scaling in (II. ip 

x := fix, t := Xi, u := vu, with 



du hL§2L &u v du d 2 u v d 2 u 

dt A dt ' dxi fj, dxi ' dx'f "J? dx'f ' 



(2.3) 

and substituting those expressions in (11. ip yields 



To keep this equation invariant, the following must be fulfilled: 



(2.4) t = so that 



// := X 13 =^> v := X fi ™ and u(x, t) := A ^ u(x, t) = X 8 n u{-, |) 
Consequently, 



u(x,t) :=t 4 ' 3 n 1 /(^), 
where t = X and /(^) = u(-%?, 1). Owing to (12. 4p . we obtain 

na + 4(3 = 1, 

which links the parameters a and (3. Hence, substituting 

(2.5) u(x,t) :=t- a f(y), with y := f , /? = 

into (11.11) and rearranging terms, we find that the function / solves a quasilinear elliptic 
equation of the form 

(2.6) V-(|/| n VA/) = Qt/ + /fyV-/. 

Finally, thanks to the above relation between a and /3, we find a nonlinear eigenvalue 
problem of the form 



l—na 
4 ' 



(2.7) _v.(|/|»VA/) + ±^yV-/ + a/ = 0, / 6 C (R N ) , 



where we add to the equation (12. 6p a natural assumption that / must be compactly 
supported (and, of course, sufficiently smooth at the interface, which is an accompanying 
question to be discussed as well). 

Thus, for such degenerate elliptic equations, the functional setting in ( 12 .7p assumes 
that we are looking for (weak) compactly supported solutions f(y) as certain "nonlinear 
eigenf unctions" that hopefully occur for special values of nonlinear eigenvalues {o;fc}fc>o- 
Our goal is to justify that, labelling the eigenfunctions via a multiindex a, 



(2.8) ( 12.71) possesses a countable set of eigenf unction/ value pairs {f a , ack}\a\=k>o- 



Concerning the well-known properties of finite propagation for TFEs, we refer to papers 
[16J, where a large amount of earlier references is available; see also [22] for more 



recent results and references in this elliptic area. However, there are still a little of 
entirely rigorous results, especially those that are attributed to the Cauchy problem for 
TFEs. In the linear case n = 0, the condition / e Co(®L N ), is naturally replaced by the 
requirement that the eigenf unctions ipp{y) exhibit typical exponential decay at infinity, 
a property that is reinforced by introducing appropriate weighted L 2 -spaces. Actually, 
using the homotopy limit n — > + , we will be obliged for small n > 0, instead of Co-setting 
in (12. 7p . we will use the following weighted L 2 -space: 

(2.9) / G L 2 p (R N ), where p(y) = e a|i/|4/3 , a > small. 

Note that, in the case of the Cauchy problem with conservation of mass making use of 
the self-similar solutions (12.51) . we have that 

M(t) := / u(x,t)dx = t~ a J f(^)dx = t- a+f3N J f(y)dy, 

R N R N ' R N 

where the actual integration is performed over the support supp / of the nonlinear eigen- 
function. Then, as is well known, if J f ^ 0, the exponents are calculated giving the first 
explicit nonlinear eigenvalue: 

(2.10) -a + f3N = a ^ n )=N an d /3 ( n ) = j-L-. 



3. Hermitian spectral theory of the linear rescaled operators 

In this section, we establish the spectrum cx(B) of the linear operator B obtained from 
the rescaling of the linear counterpart of (11.11) . the bi-harmonic equation (11.31) . which will 
be essentially used in what follows. 

3.1. How the operator B appears: a linear eigenvalue problem. Let u(x,t) be 

the unique solution of the CP for the linear parabolic bi-harmonic equation (11.31) with the 
initial data (the space as in (12. 9ft to be more properly introduced shortly) 

u e L 2 p (R N ), 

given by the convolution Poisson-type integral 

(3.1) u(x,t) = b(t) * u = t~T f F((x- z)t-i)u {z)dz. 

R N 

Here, by scaling invariance of the problem, in a similar way as was done in the previous 
section for (11.11) . the unique the fundamental solution of the operator + A 2 has the 
self-similar structure 

(3.2) b(x,t)=t-T F (y), y:=-^ (x 6 R N ). 

Substituting b(x, t) into (11.31) yields that the rescaled fundamental kernel F in (13.21) solves 
the linear elliptic problem 

(3.3) BF = — A 2 F +\y ■ W y F + ^ F = in R N , jF(y)dy=l. 




B is a non-symmetric linear operator, which is bounded from H^(R N ) to L 2 p (R N ) with 
the exponential weight as in (I2.9p . Here, a G (0, 2d) is any positive constant, depending 
on the parameter d > 0, which characterises the exponential decay of the kernel F(y): 

(3.4) \F(y)\ < Z)e-* |4/3 in R N (D > 0, d = 3 • 2~t). 

Later on, by F we denote the oscillatory rescaled kernel as the only solution of (13. 3p . 
which has exponential decay, oscillates as \y\ — > oo, and satisfies the standard pointwise 
estimate (13.41) . 

Thus, we need to solve the corresponding linear eigenvalue problem: 



(3.5) B^ = XiP in R N , ip e L 2 (R N ). 



One can see that the nonlinear one (12. TP formally reduces to ( 1 3 . 5 f) at n = with the 
following shifting of the corresponding eigenvalues: 



A = -a + f 

In fact, this is the main reason to calling (12. Tl) a nonlinear eigenvalue problem, and, 
crucially, the discreteness of the real spectrum of the linear one (13.51) will be shown to be 
inherited by the nonlinear problem, but a long way is needed to justify such an issue. 

3.2. Functional setting and semigroup expansion. Thus, we solve ( 13 .51) and calcu- 
late the spectrum of cr(B) in the weighted space L 2 (R N ). We then need the following 
Hilbert space: 



H}(R N ) CL 2 p (R N ) cl 2 * v 

r4nn>iv\ 



The Hilbert space H p (R ) has the following inner product: 

l 



(v,w) := J p(y) J2 D k v(y)D k w(y)dy, 



k=0 



where D k v stands for the vector {D@v , \/3\ = k}, and the norm 



12 ' " / P(v) E \D k v(y)\ 2 dy. 

Next, introducing the rescaled variables 

(3.6) u(x,t) =t~T W (y,T), y := -fc, r = \nt:R + ^R, 
we find that the rescaled solution w(y, r) satisfies the evolution equation 

(3.7) w T = Bw , 

since, substituting the representation of u(x,t) (I3.6P into (11.31) yields 

-Alw + {y.V y w + %w = t^. 
Thus, to keep this invariant, the following should be satisfied: 

t ^ = 1 =^> r — hit, i.e., as defined in (13.61) . 
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Hence, w(y, r) is the solution of the Cauchy problem for the equation (13.71) and with the 
following initial condition at r = 0, i.e., at t — 1: 

(3.8) w (y) = u(y, 1) = 6(1) * u = F * u . 

Thus, the linear operator — B is also a rescaled version of the standard parabolic 
one Jj; + A 2 . Therefore, the corresponding semigroup e Bl ~ admits an explicit integral 
representation. This helps to establish some properties of the operator B and describes 
other evolution features of the linear flow. From (13.11) we find the following explicit 
representation of the semigroup: 

(3.9) w(y,r)= J F(y — ze~i) Uo(z) dz = e Br w , where x = t^y } r = Int. 

R N 

Subsequently, consider Taylor's power series of the analytic kernefl 

(3.10) F(y-ze^)=Ze-^^DPF(y)z^Ze-^^mMy)z^ 

03) 03) 

for any y G Mr, where 

r /3 _/3i ~Pn 

*> ■ — A \ ' ' ' 

and ip £ are the normalized eigenf unctions for the operator B. The series in (13.101) converges 
uniformly on compact subsets in z G M N . Indeed, denoting \/3\ = I and estimating the 
coefficients 

| £ tfL D f 3 F(y)z^---z^\ < 6,1*1', 
101=' 

by Stirling's formula we have that, for I ^> 1, 

(3.11) k = fsup yeRNm=l \D?F{y)\ « f l' 1 ^ « r^V = e - /ln3/ / 4+nnc . 
Note that, the series 

Y,bi\z\ l 

has the radius of convergence R = oo. 

Thus, we obtain the following representation of the solution: 

(3.12) w(y 1 r) = Ze-^ T M l3 (u )My)^ where Xp =: 

and {^/j} are the eigenvalues and eigenfunctions of the operator B, respectively, and 

Mpiuo) :=^Jz^---z^u (z)dz 

are the corresponding momenta of the initial datum wo defined by ( 13.81) . 



1 We hope that returning here to the multiindex (3 instead of a in (|2.8| will not lead to a confusion 
with the exponent (3 in self-similar scaling (|2.5[) . 

10 



3.3. Main spectral properties of the pair {B, B*}. Thus, the following holds [T2] : 

Theorem 3.1. (i) The spectrum o/B comprises real eigenvalues only with the form 

(3.13) a(B):={A /3 = -f ,|/3|=0,1,2,...}. 
Eigenvalues \@ have finite multiplicity with eigenf unctions, 

(3.14) M y) -= { -^^F(y) = ^(4) Sl ■ ■ ■ (^Y"F(y). 

(ii) The subset of eigenf unctions $ = {ipp} is complete in L 2 (R N ) and in L 2 (R N ). 

(iii) For any A ^ cr(B), the resolvent (B — A/) -1 is a compact operator in L 2 (R N ). 

Subsequently, it was also shown in [12] that the adjoint (in the dual metric of L 2 (R N )) 
operator of B given by 

(3.15) B*:=-A 2 -iyV, 

in the weighted space L 2 ,(R N ), with the exponentially decaying weight function 

(3.16) p *(y) = ^ = e - a \yr >0 , 
is a bounded linear operator, 

B* : Hp*(R N ) -> L 2 *(R N ), so (Bv,w) = (v,B*w), v G H$(R N ), w e H 4 p ,{R N ). 

Moreover, the following theorem establishes the spectral properties of the adjoint operator 
which will be very similar to those shown in Theorem l3.1l for the operator B. 

Theorem 3.2. (i) The spectrum o/B* consists of eigenvalues of finite multiplicity, 

(3.17) cr(B*) = <r(B) := {A^ = -f , \(3\ = 0, 1, 2, ...}, 

and the eigenf unctions ip%(y) are polynomials of order \/3\. 

(ii) The subset of eigenf unctions <3>* = {ipp} is complete in L 2 , (R N ) . 

(iii) For any A ^ o~(B*) ; the resolvent (B* — A/) -1 is a compact operator in L 2 *(R N ). 

It should be pointed out that, since i/; = F and = 1, we have 

M,,Vo>= J *Pody= J F(y)dy=l. 

However, thanks to (I3.14p . we have that 

J ^ = (^,r ) = for any \(3\ ^ 0. 

This expresses the orthogonality property to the adjoint eigenf unctions in terms of the 
dual inner product. 

Note that [12] . for the eigenfunctions {ipp} of B denoted by (I3.14p . the corresponding 
adjoint eigenfunctions are generalized Hermite polynomials given by 

r [1/31/4] , 

(3-18) r p {y):= ^ £ jf A 2 V • 
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Hence, the orthonormality condition holds 

(3.19) (ipp, V> 7 ) = Spy for any ft, 7, 

where (•, •) is the duality product in L 2 {R N ) and <5g 7 is Kronecker's delta. Also, operators 
B and B* have zero Morse index (no eigenvalues with positive real parts are available). 
Key spectral results can be extended [12J to 2mth-order linear poly-harmonic flows 

(3.20) u t = -(-A) m u in R N x R + , 

where the elliptic equation for the rescaled kernel F(y) takes the form 

(3.21) BF = -(-A y rF+^yV y F+£ l F = in R N , jF(y)dy=l. 

R N 

In particular, for m = 1, we find the Hermite operator and the Gaussian kernel (see [8] 
for further information) 

BF = AF + \y • VF +f F=0 F(y) = J ^ m e-^, 

whose name is connected with fundamental works of Charles Hermite on orthogonal poly- 
nomials {Hp} about 1870. These classic Hermite polynomials are obtained by differenti- 
ating the Gaussian: up to normalization constants, 

(3.22) D^e- 1 ^ = Hp{y) e"^ for any (3. 

Note that, for N = 1, such operators and polynomial eigenfunctions in ID were stud- 
ied earlier by Jacques C.F. Sturm in 1836; on this history and Sturm's main original 
calculations, see (THl Ch. 1]. 

The generating formula f 1 3 . 2 2 j) for (generalized) Hermite polynomials is not available if 
m > 2, so that (13.181) are obtained via a different procedure, [12J. 



4. Similarity profiles for the Cauchy problem via n-BRANCHiNG 

4.1. Derivation of the branching equation. In general, construction of oscillatory 
similarity solutions of the Cauchy problem for the TFE-4 (II .ip is a difficult nonlinear 
problem, which is harder than for the corresponding FBP one. On the other hand, for 
n — 0, such similarity profiles exist and are given by eigenfunctions {ipp}- In particular, 
the first mass-preserving profile is just the rescaled kernel F(y), so it is unique, as was 
shown in Section [31 Hence, somehow, there can be expected a possibility to visualize such 
an oscillatory first "nonlinear eigenfunction" f(y) of changing sign, which satisfies the 
nonlinear eigenvalue problem (12.71) . at least, for sufficiently small n > 0. This assumes 
using the n-branching approach that "connects" / with the rescaled fundamental profile 
F satisfying the corresponding linear equation (13. 3p . with all the necessary properties of 
F presented in Section [3J 

Thus, we plan to describe the behaviour of the similarity profiles {fp}, as nonlinear 
eigenfunctions of (12. 7p for the TFE performing a "homotopic" approach when n { 0. 
Homotopic approaches are well-known in the theory of vector fields, degree, and nonlinear 
operator theory (see [TTJ [32] for details). However, we shall be less precise in order to 
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apply that approach, and here a "homotopic path" just declares existence of a continuous 
connection (a curve) of solutions / G Co that ends up at n = + at the linear eigenfunction 
ipo(y) = F(y) or further eigenfunctions ij)p(y) ~ D^F(y), as (13.141) claims. 

Using classical branching theory in the case of finite regularity of nonlinear operators 
involved, we formally show that the necessary orthogonality condition holds deriving the 
corresponding Lyapunov- Schmidt branching equation. We will try to be as much rigorous 
as possible in supporting of delivering the nonlinear eigenvalues {a k }. Further extensions 
of solutions for non-small n > require a novel essentially non-local technique of such 
nonlinear analysis, which remains an open problem. 

Those critical eigenvalues {a k } are obtained according to non-self-adjoint spectral the- 
ory from Section [3J We then use the explicit expressions for the eigenvalues and eigen- 
functions of the linear eigenvalue problem (13. 5p given in Theorem 13.11 where we also need 
the main conclusions of the "adjoint" Theorem l3.2l Then, taking the corresponding linear 
equation from (12.71) with n = 0, we find, at least, formally, that 

„ = : C(a)f := -A 2 / + \ y ■ V/ + af = 0. 

From that equation combined with the eigenvalues expressions obtained in the previous 
section, we ascertain the following critical values for the parameter a k = a k (n), 

(4.1) n = 0: a k (0):=-X k + ^ = ^f for any k = 0, 1, 2, . . . , 
where X k are the eigenvalues defined in Theorem l3.ll so that 

a (0) = f , ttl (0) = ^i, 02(0) = . . . , a fc (0) = . . . . 
In particular, when k = 0, we have that a (0) = % and the eigenfunction satisfies 
BF = 0, so that ker C(q.q) = span {^o} (i ) o — F'), 

and, hence, since A = is a simple eigenvalue for the operator £(oq) = B, its algebraic 
multiplicity is 1. In general, we find that 

(4.2) ker (B + | J) = span {ipp, \(3\ = k}, for any k = 0, 1, 2, 3, • • • , 

where the operator B + | / is Fredholm of index zero since it is a compact perturbation of 
the identity of linear type with respect to k. In other words, R[C(otk)] is a closed subspace 
of Lp(R. N ) and, for each a k , 

dimker(£(afc)) < oo and codim_R[£(afe)] < oo. 
Then, for small n > in (12.71) . we can use the asymptotic expansions 

(4.3) Ofc(n) := a k + Hi,k n + o(n), and 

(4.4) \f\ n = e nln|/l := 1 + nln |/| + o(n). 

As customary in bifurcation-branching theory [321 S3] , existence of an expansion such as 
( 14. 3 p will allow one to get further expansion coefficients in 

(4.5) Ofe(n) := a k + // 1(fc n + \i^ k n 2 + ii^ k r? + ... , 
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as the regularity of nonlinearities allows and suggests, though the convergence of such an 
analytic series can be questionable and is not under scrutiny here. 

Another principle question is that, for oscillatory sign changing profiles f(y), the last 
expansion (14.41) cannot be understood in the pointwise sense, but can be naturally ex- 
pected to be valid in other metrics such as weighted L 2 or Sobolev spaces as in Section [3] 
that used to be appropriate for the functional setting of the equivalent integral equation 
and for that with n = 0. Since (I4.4p is obviously pointwise violated at the nodal set 
{/ = 0} of f(y), this imposes some restrictions on the behaviour of corresponding eigen- 
functions ipp{y) (n = 0) close to their zero sets. Using well-known asymptotic and other 
related properties of the radial analytic rescaled kernel F(y) of the fundamental solutions 
(13. 2p . the generating formula of eigenfunctions (I3.14p confirms that the nodal set of ana- 
lytic eigenfunctions {ipp = 0} consists of isolated zero surfaces, which are "transversal" , 
at least in the a.e. sense, with the only accumulation point at y = oo. Overall, under 
such conditions, this indicates that 

(4.6) expansion f!4.4j) contains not more than "logarithmic" singularities a.e., 

which well suited the integral compact operators involved into a branching analysis, 
though we are far away to claim this as any rigorous issue. Moreover, when n > is 
not small enough, such an analogy and statements like (14.61) become not that clear, and 
global extensions of continuous n-branches induced by some compact integral operators, 
i.e., nonexistence of turning (saddle-node) points in n, require, as usual, some unknown 
monotonicity-like results. 

Now, we assume the expansion (14.41) away from possible zero surfaces of f{y), which, by 
transversality, can be localized in arbitrarily small neighbourhoods. Indeed, it is clear that 
when \f\ > 5 > 0, for any S > 0, there is no problem in approximating \f\ n by (14. 3p . i.e., 
|/| n = 1 + 0(n) as n — ► + . However, when |/| < 5, for any 5 > sufficiently small, the 
proof of such an approximation in weak topology (as suffices for dealing with equivalent 
integral equations) is far from clear unless the zeros of the /'s are also transversal a.e. 
with a standard accumulating property at the only interface zero surface. The latter issues 
have been studied and described in [TJ] in the radial setting. Hence, we can suppose that 
such nonlinear eigenfunctions f(y) are oscillatory and infinitely sign changing close to the 
interface surface. Therefore, if we assume that their zero surface are transversal a.e. with 
a known geometric-like accumulation at the interface, we find that, for any n close to zero 
and any 5 = 5(n) > sufficiently small, 

n|ln|/||»l, if \f\<5(n), 

and, hence, on such subsets, f(y) must be exponentially small: 

|ln|/||»i =► ln|/|«-i =► |/|«e-i 

Recall that this happens in also exponentially small neighbourhoods of the transversal 
zero surfaces. 

Overall, using the periodic structure of the oscillatory component at the interface [H] 
(we must admit that such delicate properties of oscillatory structures of solutions are 
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known for the ID and radial cases only, though we expect that these phenomena are 
generic), we can control the singular coefficients in (14.31) . and, in particular, to see that 

(4.7) H/| G LtoJR"). 

However, for most general geometric configurations of nonlinear eigenf unctions f(y), we 
do not have a proper proof of ( 14. 7ft or similar estimates, so our further analysis is still 
essentially formal. It is worth recalling again that our computations below are to be 
understood as those dealing with the equivalent integral equations and operators, so, in 
particular, we can use the powerful facts on compactness of the resolvent (B — A/) -1 
and of the adjoint one (B* — A/) -1 in the corresponding weighted L 2 -spaces. Note that, 
in such an equivalent integral representation, the singular term in (14.41) satisfying (14.71) 
makes no principal difficulty, so the expansion (14.41) makes rather usual sense for applying 
standard nonlinear operator theory. 

Thus, under natural assumptions, substituting (14. 3 j) into (12.71) . for any k = 0, 1, 2, 3, • • • , 
we find that, omitting o(n) terms when necessary, 

-V • [(1 + n In |/|) VA/] + ^Y 1 '^ 2 V ■ W + (a* + /i M n)/ = , 
and, rearranging terms, 

-A 2 / - nV • (In |/| VA/) + \y - V/ - ^±p^ y . V/ + a k f + ^ k nf = . 
Hence, we finally have 

(4.8) (B + f /)/ + n[ - V • (In |/| VA/) - ■ V/ + ^, k f] + o{n) = , 
which can be written in the following form: 

(4.9) (B + f/)/ + aA4(/) + o(n) = 0, 
with the operator 

(4.10) Af k (f) := -V ■ (In |/| VA/) - ^ y - V/ + ^, k f . 
Subsequently, as was shown in Section [31 we have that 

(4.11) ker (B + | /) = span {iffp}\p\=h for any k = 0, 1, 2, 3, • • • , 
where the operator B + 1 1 is Fredholm of index zero and 

dimker (B + \ I) = M k > 1 for any jfe = 0, 1, 2, 3, • • • , 

where M k stands for the length of the vector {D@v, \/3\ = k}, so that M k > 1 for k > 1. 
Simple eigenvalue for k = . Since is a simple eigenvalue of B when k = 0, i.e., 

ker B®R[B] = L 2 p (R N ), 

the study of the case k = seems to be simpler than for other different fc's because the 
dimension of the eigenspace is M = 1. Thus, we shall describe the behaviour of solutions 
for small n > and apply the classical Lyapunov-Schmidt method to (14.91) (assuming 
as usual some extra necessary regularity hypothesis to be clarified later on), in order to 
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accomplish the branching approach as n [ 0, in two steps, when k = and k is different 
from 0. 

Thus, owing to Section [3j we have already known that is a simple eigenvalue of B, 
i.e., ker B = spanj-^o} is one-dimensional. Hence, denoting by Yq the complementary 
invariant subspace, orthogonal to i/jq, we set 

/ = + V , 

where Vq G Yq. According to the already well known spectral properties of the operator 
B, we define Pq and P\ such that Pq + P\ = /, to be the projections onto ker B and Yq 
respectively. Finally, setting 

(4.12) V Q :=n$i >0 + o(n), 

substituting the expression for / into (14.91) and passing to the limit as n — > + leads to a 
linear inhomogeneous equation for $1,0, 

(4.13) B$i, = -JVo(Vo), 

since B?/>o = 0. Moreover, by Fredholm theory, Vq G Yq exists if and only if the right-hand 
side is orthogonal to the one dimensional kernel of the adjoint operator B* with ip^ = 1, 
because of (I3.18p . Hence, in the topology of the dual space L 2 , this requires the standard 
orthogonality condition: 

(4.14) <Aro(^o),l> = 0. 

Then, (14. 13|) has a unique solution $ 10 G Yq determining by (14.121) a bifurcation branch for 
small n > 0. In fact, the algebraic equation (14.141) yields the following explicit expression 
for the coefficient yU 10 of the expansion (14.3j) of the first eigenvalue a (n): 

(4-15) ^,0 := <v-(in|.o|VA^,v, ^) = (y . (ln WAi;o) + N y . 

Multiple eigenvalues for k > 1 . For any k > 1, we know that 

dim ker (B + \ I) = M k > 1. 
Hence, we have to use the full eigenspace expansion 

(4.16) /= J2 C/#/3 + V k , 

\p\=h 

for every k > 1. Currently, for convenience, we denote {tp/3}\i3\ = k = {^1, ipM k } the 
natural basis of the M^-dimensional eigenspace ker (B + 1 1) and set ipk = Y^\/3\=k c P' l l ) P- 
Moreover, Vk G Yk and Vk = Tl l \p\>k c P' i l J P^ where Y k is the complementary invariant 
subspace of ker (B + | /) . Furthermore, in the same way, as we did for the case k = 0, 
we define the Po,fc and P± : k, for every k > 1, to be the projections of ker (B + | /) and 
respectively. We also expand V k as 

(4.17) V k :=n$i )fc + o(n). 
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Subsequently, substituting (14.161) into (14.91) and passing to the limit as n [ + , we obtain 
the following equation: 

(4-18) (B + J J)*i, fc = -MT,\p\= k ^p), 

under the natural "normalizing" constraint 

(4.19) E C/3 = 1 ic p > 0). 

|/3|=fc 

Therefore, applying the Fredholm alternative, V). G Y k exists if and only if the term on 
the right-hand side of (I4.18P is orthogonal to ker (B + \ 1) . Multiplying the right-hand 
side of ( 14.18P by for every \/3\ = k, in the topology of the dual space L 2 , we obtain an 
algebraic system of + 1 equations and the same number of unknowns, {cp, |/3| = k} 
and yUi jfe : 

(4.20) (M k (.E\0\=k c M,ty) = O for all \(3\ = k, 

which is indeed the Lyapunov-Schmidt branching equation [43] . In general, such algebraic 
system are assumed to allow us to obtain the branching parameters and hence establish 
the number of different solutions induced on the given M^- dimensional eigenspace as the 
kernel of the operator involved. 

However, we must admit and urge that the algebraic system (14.201) is a truly difficult 
issue. One of the main features of it is as follows: 



(4.21) fECffll) is not variational 



In other words, one cannot use for (14.201) the classic category-genus theory of calculus of 
variation [U[32], to claim that the category of the kernel (equal to M^) is the least number 
of different critical points and hence of different solutions. 

To see (14.211) . it suffices to note that, due to (13.141) and (I3.18p . the generalized Hermite 
polynomials have nothing common in the algebraic sense with the eigenf unctions ipp 
in the L 2 -scalar products in (I4.20p . 

4.2. A digression to Hermite classic self-adjoint theory. It is worth mentioning 
that, for classic second-order Hermite operator 

(4.22) B = A + |y-V + fJ (then, in the L 2 -metric, B* = A - \ y ■ V) , 

(14.211) is not the case. Indeed, by classic self-adjoint theory [H p. 48], these eigenfunctions 
are related to each other by 

(4.23) My)=D p F(y) = Hp(y)F(y), where F(y) = (4tt)-t e ~— 

is the Gaussian kernel and Hp(y) are standard Hermite polynomials, which also define 
the adjoint eigenfunctions: 

(4-24) PM = b f* H M = 4faMv)i 
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where bp are normalization constants. One knows that this is a result of the symmetry of 
the operator (I4.22p in the weighted metric of L 2 (W N ), where 



P(l/) = e* =► B = IV-(pV) + f /, so (B)* ? = B. 

In view of the relations (14.231) and (I4.24p of the bi-orthonormal bases {iip} and 
the corresponding algebraic systems such as (14.201) can be variational. Moreover, even 
the original nonlinear elliptic equation similar to (12. 7p . where the 4th-order operator is 
replaced by a natural 2nd-order one of the porous medium type: 

-V(|/PVA/) V(|/|"V/), 

then becomes variational itself. Thus, in this case, both branching (local phenomena) 
and global extensions of n-bifurcation branches can be performed on the basis of pow- 
erful Lusternik-Schnirel'man category variational theory from 1920s [221 § 56], so that 
existence and multiplicity (at least, not less than in the linear case n = 0) of solutions are 
guaranteed. 

4.3. Computations for branching of dipole solutions in 2D. To avoid excessive 
computations and as a self-contained example, we now ascertain some expressions for 
those coefficients in the case when = 1, N — 2, and Mi = 2, so that, in our notations, 
{ipp}\/3\=x = {ipi, ip 2 }- Consequently, in this case, we obtain the following algebraic system: 
the expansion coefficients of ipi = ci^i + c 2 ?/> 2 satisfy 

[ cSl, hi) - 3f {Vx, V ■ Wl> + C lf i 1A + cSh h 2 ) - Sfl y . v4) = 0, 
(4.25) I Cl (4*, hx) - ssi (4*, y . + c 2 (r 2 , h 2 ) - Sfi (4*, y ■ V4> + ca/ii,! = 0, 
I Ci + c 2 = 1, 

where 

h x := -V • Mci4>i + ^2)VA^], h 2 : = -V • pn(cA + c 2 ^ 2 )VA^ 2 ], 

and, ci, c 2 , and /i^i are the coefficients that we want to calculate, ax is regarded as the 
value of the parameter a denoted by (14. ip and dependent on the eigenvalue Ai, for which 
ipx,2 are the associated eigenfunctions, and ^* 2 the corresponding adjoint eigenfunctions. 
Hence, substituting the expression c 2 = 1 — cx from third equation into the other two, we 
have the following nonlinear algebraic system 

= Nx(cx, - cxf [(fay ■ Wi) - (i>*x, y ■ V^>] , 
= N 2 ( Cl , M - Cl f [(4*, y • V^i) - (4% 1/ • W 2 >] + 



(4.26) 
where 



iVi(ci, /ii,i) := Ci<^J, M + /l2> - f (V»i,2/ • V^ 2 ) - Cl (^, /i 2 > + CxtMx,!, 

N 2 {ci,m tl ) := cx(^ 2 ,hx) + (ip2,h 2 ) - f $ly ■ V4> 2 ) - cx(4>lh 2 ) - cxVx,i 

represent the nonlinear parts of the algebraic system, with ho and hx depending on c\. 

Subsequently, to guarantee existence of solutions of the system ( I4.25p . we apply the 
Brouwer Fixed Point Theorem to (I4.26P by supposing that the values ci and fix,x are the 
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unknowns, in a disc sufficiently big D R (ci, /2x,i) centered in a possible nondegenerate zero 
Thus, we write the system ( I4.26P in the matrix form 

°V ( -f ft,!/'VW-ferV«] \(c x \ (N x {c x ^ lx 



.1, 



are the possible solutions of (14.261) . where 9Jt is the matrix corresponding to the linear 



Hence, we have that the zeros of the operator 

^(ci,/ii,i) :=VH' 

Hi 

part of the system, while 

(iVi(ci, fi 1A ), N 2 (c 1 , /ii,i)) T , 
corresponds to the nonlinear part. The application 7i : A x [0, 1] — > R, defined by 

H^m^t) ■.= an +t\ 



A*i,i/ V^V" 2 (ci,/ii, : 

provides us with a homotopy transformation from the function JF(ci, //x,x) = H(ci, 1) 
to its linearization 



(4.27) W(ci,/xi,i,0) :=0tt 



ci 



Thus, the system (I4.26P possesses a nontrivial solution if (14.271) has a nondegenerate 
zero, in other words, if the next condition is satisfied 

(4.28) fe r VW-(ft r V«^0. 

Note that, if the substitution would have been c\ — 1 — C2, the condition might also be 

Then, under condition (14.281) . the system (I4.26P can be written in the form 

(4.29) ( 01 ~ £ ! \=-M^(^^~ h 

which can be interpreted as a fixed point equation. Moreover, applying Brouwer's Fixed 
Point Theorem, we have that 

Ind((cx, /ix,!), H(; ; 0)) = Qc^fr.oW, '» °)) 

= Deg(W(.,.,0),iPj i (ci,/ii,i)) 
= Deg(J r (ci, y Ui i i),D /? (ci,/ii i i)), 

where Qc fl (ci,/ii,i)(^(', ■, 0)) defines the number of rotations of the function Tl(-,-,0) 
around the curve Cr(ci, /tx,x) and Deg(7i(-, •, 0), -Dr(cx, /ix,i)) denotes the topological de- 
gree of 7^(-, -,0) in Dr(ci, /ti,i). Owing to classical topological methods, both are equal. 
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Thus, once we have proved the existence of solutions, we achieve some expressions for 
the coefficients required: 



Pl.i = + fc, h - h 2 ) - ^ (fe + y-Vfa-y. W> 2 )) 



+ v 2 , M + fW+^>r Wi>, 



Cl 



c 2 . 



The expressions for the coefficients in a general case might be accomplished after some 
tedious calculations, otherwise similar to those performed above. Note that, in general, 
those nonlinear finite-dimensional algebraic problems are rather complicated, and the 
problem of an optimal estimate of the number of different solutions remains open. More- 
over, reliable multiplicity results are very difficult to obtain. We expect that this number 
should be somehow related (and even sometimes coincides) with the dimension of the 
corresponding eigenspace of the linear operators B + 1 1, for any k — 0, 1, 2, ... . This is a 
conjecture only that may be too illusive; see further supportive analysis presented below. 

However, we devote the remaining of this section to a possible answer to that conjecture, 
which is not totally complete though, since we are imposing some conditions. 

Thus, in order to detect the number of solutions of the nonlinear algebraic system 
(14.251) . we proceed to reduce this system to a single equation for one of the unknowns. 
As a first step, integrating by parts in the terms in which h± and h 2 are involved and 
rearranging terms in the first two equations of the system (14.251) . we arrive at 



/ V$l ■ hi(ci^i + c 2 ^ 2 )VA(ci^i+c 2 ^ 2 ) 

R N 

-cif J kv ■ Wi + Cl /i M - c 2 f / fay ■ V^ 2 = 0, 

R N R N 

J • ln(ci^i + c 2 ^ 2 )VA(c 2 ^i + c 2 ^ 2 ) 

R N 

-cif / kv ■ Wi + c 2fJll>1 -c 2 f J &y ■ V^ 2 = 0. 



By the third equation, we have that c\ = 1 — c 2 , and hence, setting c\ipi + c 2 ip 2 = 
■01+C02— i>i) c 2 and substituting these into those new expressions for the first two equations 
of the system, we find that 



(4.30) 



' / Vipl ■ ln(^i+(^ 2 - ^i)c 2 )VA(^ + (^ 2 - k)c 2 ) + - c 2 /i lyl 

N 

-f / fay ■ Vk + c 2 f j fay ■ (V^i - V^ 2 ) = 0, 

R N R N 

J Vfe ■ Hi>l + (k - ^l)c 2 )VA(^! + $ 2 - k)c 2 ) + C2/il,l 

l N 

-f / fay ■ W>i + c 2 ^ / fey ■ (Vk - W> 2 ) = 0. 
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Subsequently, adding both equations, we have that 

= - f (V^ + V^|) • ln($i + {fa - ^i)c 2 )VA(^ + (^j - i) X )c 2 ) 

R N 

+f I + ^2)2/ • W>i - c 2 f / (fe + • (V^ 2 - v&). 

R'v R^ 

Thus, substituting it into the second equation of (I4.30p . we obtain the following equation 
with the single unknown c 2 : 

(4.31) 

-c|f / (fa + k)v ■ (W2 - V^i) + c 2 ^( / (fa + 2fc)y ■ V^i - / 4*2/ • W> 2 ) 

R^ R^ R^ 

-f / • v^i + / vr 2 ■ M$i + ($2 - ^i)c 2 )va(^ + $ 2 - ^)c 2 ) 

-c 2 / (V$f + Vfe) • 1b(^i + $2 - ^i)c 2 )VA(4 + (ij 2 - ^)c 2 ) = 0, 
which can be written in the following way: 

(4.32) c\A + c 2 B + C + u(c 2 ) = £(c 2 ) + u{c 2 ) = 0. 

Here, uj(c 2 ) can be considered as perturbation of the quadratic form $(c 2 ) with the coef- 
ficients defined by 

A:=-f /(^ + ^I)y-(W 2 -Wi), 

R N 

B:=f(J (fe + W 2 )y ■ V^i - / fey ■ W> 2 )„ C := -f / • V^, 

R^ RJV R^ 

cu(c 2 ) := / W> 2 * ■ ln(4 + (A - V>i)c 2 )VA(^ + (^ - ^)c 2 ) 

R^ 

-ca / (Wj + Vfe) ■ ln(^ + $ 2 - ^ 1 )c 2 )VA(^ 1 + (^ - ^)c 2 ). 

R^ 

Since, due to the normalizing constraint (I4.19p . c 2 G [0, 1], solving the quadratic equa- 
tion #(c 2 ) yields: 

(i) c 2 = £(0) = C; 

(ii) c 2 = 1 5(1) = A + 5 + C; and 

(iii) differentiating J with respect to c 2 , we obtain that $'(c 2 ) = 2c 2 A + B. Then, the 
critical point of the function £ is c\ = — £ and its image is 5(c 2 ) = — jt + C. 

Consequently, the conditions that must be imposed in order to have more than one 
solution (we already know the existence of at least one solution) are as follows: 

(a) C(A + B + C) > 0; 

(b) C(- ^ + C) < 0; and 

(c) < < 1. 

Note that, for ~ta + C = 0, we have just a single solution. 
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Hence, considering the equation again in the form 



d(c 2 )+oo(c 2 ) = 0, 

where 00(02) is a perturbation of the quadratic form #(02), and bearing in mind that the 
objective is to detect the number of solutions of the system (I4.25p . we need to control 
somehow this perturbation. Under the conditions (a), (b), and (c), #(02) possesses exactly 
two solutions. Therefore, controlling the possible oscillations of the perturbation 00(02) in 
such a way that 

||w(c2)|| LOO < $(c* 2 ), 

we can assure that the number of solutions for (I4.25P is exactly two. This is the dimension 
of the kernel of the operator B + 1 I (as we expected in our more general conjecture). 

The above particular example shows how difficult are the questions on existence and 
multiplicity of solutions for such non-variational branching problems. Recall that the 
actual values of the coefficients A, B, C, and others, which the number of solutions 
crucially depend on, is difficult even estimate numerically in view of a complicated nature 
of the eigenfunctions (13.141) involved, to say nothing of the nonlinear perturbation 00(02)- 



4.4. Branching computations for \/3\ = 2. Overall, the above analysis provides us 
with some expressions for the solutions for the self-similar equation (12. 7p depending on 
the value of k. Actually, we can achieve those expressions for every critical value ak, but 
again the calculus gets rather difficult. For the sake of completeness, we now analyze 
the case \@\ = 2 and M 2 = 3, so that {^g}^^ = {1^1, if>2, ^3} stands for a basis of the 
eigenspace ker (B + with k = 2 (\ k = — |). 

Thus, in this case, performing in a similar way as was done for (14.251) with ip 2 = 
Ciipi + c 2 ip2 + £3^3, we arrive at the following algebraic system: 

' CS*1, hi) + C 2 {fa h 2 ) + C 3 {fa h 3 ) - (fa y . Vfa - (fa y . V ^ 2 ) 

-SfL(fay-Vik) + OL^ = 0, 



(4.33) < 



ci(i%, h) + 02(^*2, h 2 ) + c 2 (r 2 , h) - Sf (ft, y ■ Wi) - <Va> V ' W<2> 
Cl (fa h) + c 2 (fa h 2 ) + c 2 (fa h 3 ) - (fa y ■ Vfa - (fa y ■ V^ 2 ) 

- a f a <feyV^ 3 ) + c 3/ x M = o, 

Cl + c 2 + c 3 = 1, 



where 



hi := -V ■ pn(ci^i + c 2 ip2 + CswVA^i], 

/i 2 := -V • [ln(ci^i + c 2 ^ 2 + c 3 4)VA^ 2 ], 

/13 : = -V • [ln(ci^i + c 2 ^ 2 + c 3 4)VA4], 
and Ci, c 2 , c 3 , and /xi )2 are the unknowns to be evaluated. Also, a 2 is regarded as the 
value of the parameter a denoted by (14.11) and is dependent on the eigenvalue A 2 with 
i>2, i>3 representing the associated eigenfunctions and t/j*, i/;^, ^3 the corresponding ad- 
joint eigenfunctions. 
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Subsequently, substituting c 3 = 1 — c x — c 2 into the first three equations and performing 
an argument based upon the Brouwer Fixed Point Theorem and the topological degree 
as the one done above for the case |/3| = 1, we ascertain the existence of a nondegenerate 
solution of the algebraic system if the following condition is satisfied: 

(4.34) (^^y•V(4-^ 1 )>(^2 J ^/•V(^3-^2))-(^I,y•V(4-A))<^^l/•V(^3-^l)>^0. 

Note that, by similar substitutions, other conditions might be obtained. 

Furthermore, once we know the existence of at least one solution, we proceed now with 
a possible way of computing the number of solutions of the nonlinear algebraic system 
(14.331) . Obviously, since the dimension of the eigenspace is bigger than that in the case 
when |/3 1 = 1, the difficulty to obtain multiplicity results increases. 

Firstly, integrating by parts in the nonlinear terms, in which hi, h 2 and h 3 are involved, 
and rearranging terms in the first three equations gives 

/ VV>* • ln(ci^i + c 2 ^ 2 + c 3 f/> 3 )VA(ci^i + c 2 ^ 2 + c 3 ^ 3 )-ci^ / 4>*y ■ Wi + ci//i, 2 



R N 



-c 2 f / kv ■ W> 2 - c 3 f / ■ = 0, 



/ • hi(ci^i + c 2 ^ 2 + c 3 4)VA(c 2 fi + c 2 ^ 2 + c 3 ^ 3 )-ci^ / fay ■ V^i + c 2 /ii, 2 

R N R N 

-c 2 f / fy ■ V0 2 - c 3 f / fay ■ V4 = 0, 

R N R N 

f • ln(ci^i + c 2 2 + C3tA 3 )VA(c 2 ^i + c 2 ^ 2 + c# 3 )-Ci^ / ifey ■ V^i + c 3y ui i2 

R N R N 

-c 2 f / kv ■ V^ 2 - c 3 f / r 3 y ■ v4 = o. 

R N R N 

By the fourth equation, we have that c\ = 1 — c 2 — c 3 . Then, setting 

Cl-01 + C 2 ^2 + C 3 ljj 3 = V>i + C 2 (?/>2 - -01 ) + C 3 (?/> 3 - 

and substituting it into the expressions obtained above for the first three equations of the 
system yield 

/ V^J • ln(^i+(0 2 - ^i)c 2 + (4 - ^)c 3 )VA(^ + (0 2 - ^i)c 2 + (4 - ^i)c 3 ) 

R^ 

+/il,2 - C 2 //l, 2 - C 3 /ii i2 - X / ^ ' 

R^ 

+f / 01*2/ ■ ((W>i - V^ 2 )c 2 + (V^a - V4)c 3 ) = 0, 

R N 
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/ V4 ■ hi(4+(4 - 4)c 2 + (A - *)c 3 )VA(^ + (4 - 4)c 2 + (4 - 4)c 3 ) 

R N 

(4 35) +c 2 /ii, 2 - x / 4?/ • Wi 

+f / 4*2/ • ((v4 - v4)c 2 + (v4 - v4)c 3 ) = o, 

R N 

J V4 ■ ln(4+(4 - 4)c 2 + (4 - ^i)c 3 )VA(^ + (4 - 4)c 2 + (4 - 4)c 3 ) 

R N 

+C3//1.2 - f / 4*2/ • Wi 

+f / 4*2/ ■ ((Wi - V4)c 2 + (v4 - v4)c 3 ) = o. 

R N 

Now, adding the first equation of (14.351) to the other two, we have that 
/ (V4 + V4) • ln(4+(4 - 4)c 2 + (4 - 4)c 3 )VA(4 + (4 - $1)02 + (4 - 4)c 3 ) 

R^ 

+/ii, 2 - c 3 /i lj2 - f / (4 + 4)y ■ V4 

R^ 

+f / (4 + 4)2/ ■ ((Wx - v4)c 2 + (V4 - v4)c 3 ) = o, 

R^ 

/ (V4 + V4) • hi(4+(4 - 4)c 2 + (A - 4)c 3 )VA(4 + (4 - 4)c 2 + (A - 4)c 3 ) 

R N 

+/il,2 - 0201,2 - f / (4 + 4)2/ • V4 

R^ 

+f / (4 + 4)2/ ■ ((v4 - v4)c 2 + (v4 - v4)c 3 ) = o. 

R^ 

Subsequently, subtracting those equations yields 

01,2= ^ [ / (V4 - V4) ■ lnM>VAM> - f / (4 - 1>*)y ■ V4 

R N R N 

+f / (4 - 4)2/ ■ ((v4 - v4)c 2 + (v4 - v4)c 3 )] , 



where \I/ = 4 + (4 — 4) c 2 + (4 — 4) c 3- Thus, substituting it into (14.351) (note that, 
from the substitution into one of the last two equations, we obtain the same equation), 
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we arrive at the following system, with c 2 and c 3 as the unknowns: 

c 3 / (W* - W 2 *+V^) • In tfVAtf - c 2 / (V^J + - V$J) • IntfVAtf 

R N R N 

+ / - V4*) • In * - f / (4* - • 

R'v R^ 

+c 2 f [ / (k - k)v ■ v(2^ - k) - J fcv ■ W>i] 

R N R N 

+c 3 f [ / (k - k)v ■ v(2^i - 4) - / • v^i] 

R N R N 

+c 2 c 3 f [ / • (V4 - W> 2 ) - / (VV> 2 * - W>*)y • (2V^ - W 2 - W> 3 )] 

R N R N 

+4f J (k -k + k)y ■ (v^i - vk) 

R N 

-4f J (k + k- k)y ■ (v^i - v4) = o, 

R^ 

c 3 / Vk ■ In WAtf -c 2 J Vk ■ In WAtf - c 3 f / 4*1/ • V^i + c 2 f / 4*1/ • 

R^ R^ 1" R N 

+c 3 f / • ((v4 - v4)c 2 + (v^i - v4)c 3 ) 

R^ 

-c 2 f / ky ■ ((v4 - v4)c 2 + (v^i - v4)c 3 ) = o. 



These can be re-written in the following form: 

Aic\ + Bic\ + dc 2 + -Dic 3 + -Eic 2 c 3 + c^i(c 2 , c 3 ) = 0, 
A 2 c 2 + B 2 c\ + C 2 c 2 + _D 2 c 3 + E 2 c 2 C3 + cj 2 (c 2 , c 3 ) = 0, 



(4.36) 
where 



and 



wi(c2, c 3 ):= c 3 / (V4 - V4 + V^) • In #VA# 

R^ 

-c 2 / (V4 + V4 - V^) • IntfVAtf 

R N 

+ I ( W 2 * - v4) • in * - f / (4 -k)y V^i 

R^ R^ 

^ 2 (c 2 ,c 3 ):= c 3 / • In $VA$ - c 2 / V^*dnWA* 



are the perturbations of the quadratic polynomials 

$i(c 2 , c 3 ) := AjC 2 + Bicl + CiC 2 + Ac 3 + EiC 2 c 3 , with i = 1, 2. 
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The coefficients of those quadratic expressions are given by 



Bi := f / (A - 4* + • (vA - v4), 

R^ 

Ci:=f[J (r 2 - k)y ■ V(2& - 4) - / ^12/ ■ Wi], 

R^ R^ 
R w R^ 

:= f [ / fe/ • (W> 3 - W> 2 ) - / (W 2 * - V^*)y ■ - - V^)], 

M'V R^ 

^:=-f Jivrt-v*), 

R^ 

S 2 :=f /^•(Wi-V^) I 

R^ 

R^ 

D 2 :=-f ffey-Vfa, 

R^ 

E 2 :=ff r 2 y ■ (V^ a - V^ 2 ) - ^ • (V& - V^ 3 )- 

R^ 

Therefore, using the conic classification to solve (14.361) . we will have the number of 
solutions through the intersection of two conies. Then, depending on the type of conic, 
we shall always obtain one to four possible solutions for our system. Hence, somehow, the 
number of solutions depends on the coefficients we have for the system and, at the same 
time, on the eigenfunctions that generate the subspace ker (B + |) . Thus, we have the 
following conditions, which will provide us with the conic section of each equation of the 
system (14.361) : 

(i) If Bf—AAiEi < 0, the equation represents an ellipse, unless the conic is degenerate, 
for example c\ + c\ + k = for some positive constant k. So, if Ai = Bj and Ei = 0, 
the equation represents a circle; 

(ii) If Bf — 4AiEi = 0, the equation represents a parabola; 

(iii) If Bf —AAiEi > 0, the equation represents a hyperbola. If we also have Ai + Ei = 
the equation represents a hyperbola (a rectangular hyperbola). 

Consequently, the zeros of the system (I4.36P and, hence, of the system (I4.33p . adding the 
"normalizing" constraint (I4.19p . are ascertained by the intersection of those two conies in 
(14.361) providing us with the number of possible n-branches between one and four. Note 
that in case those conies are two circles we only have two intersection points at most. 
Moreover, due to the dimension of the eigenspaces it looks like in this case that we have 
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four possible intersection points two of them will coincide. However, the justification for 
this is far from clear. 

Moreover, as was done for the previous case when \(3\ = 1, we need to control the 
oscillations of the perturbation functions in order to maintain the number of solutions. 
Therefore, imposing that 

\\ui(c 2 ,c 3 )\\ Loo <$i(c%,c* 3 ), with i = l,2, 

we ascertain that the number of solutions must be between one and four. This again gives 
us an idea of the difficulty of more general multiplicity results. 

4.5. Further comments on mathematical justification of existence. We return to 
the self-similar nonlinear eigenvalue problem (12.71) . associated with (jl.ip . which can be 
written in the form 

(4.37) C(a,n)f + Af(n,f) = 0, where M(n, f) := V • ((1 - |/|*)VA/) . 

As we have seen, the main difficulty in justifying the n-branching behaviour concerns 
the distribution and "transversal topology" of zero surfaces of solutions close to finite 
interface hyper-surfaces. 

Recall that, as in classic nonlinear operator theory [TTj 1321 |4~3] . our analysis above 
always assumed that we actually dealt with and performed computations for the integral 
equation: 

(4.38) / = n)Af(n, f) = Q{n, /), C(a, n) := -A 2 + ±=f* y • V + al, 

where C(a, n) is invertible in L 2 (this is directly checked via Section[3]) and, hence compact, 
for a fixed a, and / 6 Cq(M. n ) for small n > 0. This confirms that the zeros of the function 
Tin, f) are fixed points of the map Q(n,f). Note again that (I4.38P is an eigenvalue 
problem, where admissible real values of a are supposed to be defined together with its 
solvability. This makes existence/multiplicity questions for (14.381) extremely difficult. 

There are two cases of this problem. The first and simpler one occurs when the eigen- 
value a is determined a priori, e.g., in the case k = 0, where ao(O) = ^ denoted as 
Qfo(0) = a 0i an d where, for n > 0, the first nonlinear eigenvalue is given explicitly (see 
(|230])): 

«o(n) = 

Then ( 14.38)) with a = ao(n) for n > becomes a standard nonlinear integral equation 
with, however, a quite curious and hard-to-detect functional setting. Indeed, the right- 
hand side in (14.381) . where the nonlinearity is not in a fully divergent form, assumes the 
extra regularity at least such as 

(4.39) / e H 3 p . 

In view of the known good properties of the compact resolvent (C — A/) -1 , it is clear that 
the action of the inverse one C~ l is sufficient to restore the regularity, since locally in M> N 
this acts like A -2 . Therefore, it is plausible that 

(4.40) G : H 3 p - H% 
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and it is not difficult to get an a priori bound at least for small enough /'s. The accompa- 
nying analysis as y — > oo (due to the unbounded domain) assumes no novelties or special 
difficulties and is standard for such weighted L 2 and Sobolev spaces. 

Therefore, application of Schauder's Fixed Point Theorem (see e.g., jH p. 90]) to (14.381) 
is a most powerful tool to imply existence of a solution, and moreover a continuous curve 
of fixed points T n = {/, n > small}. By scaling invariance of the similarity equation, 
we are obliged to impose the normalization condition, say, 

(4.41) /(0) = 5 > sufficiently small. 

Uniqueness remains a completely open problem. However, studying the behaviour of 
the solution curve r n as n — > and applying (under suitable hypothesis) the branching 
techniques developed above, we may conclude that any such continuous curve must be 
originated at a properly scaled eigenfunction tp = F, so that such a curve is unique due 
to well-posedness of all the asymptotic expansions. 

A possibility of extension of T n for larger values of n > represents an essentially more 
difficult nonlocal open problem. Indeed, via compactness of linear operators involved in 
(14.381) . it is easy to expect that such a curve can end up at a bifurcation point only (unless 
blows up). However, nonexistence of turning saddle-node points at some > (meaning 
that the n-branch is nonexistent for some n > n*) is not that easy to rule out. Moreover, 
such turning points with thin film operators involved are actually possible, [20J. 

After establishing existence of such solutions for small n > 0, we face the next problem 
on their asymptotic properties including the fact that these are compactly supported. On 
a qualitative level, these questions were discussed in [T3] . 

In the case of higher-order nonlinear eigenfunctions of (14.381) for k > 1 including the 
dipole case k — 1, the parameter a becomes an eigenvalue that is essentially involved into 
the problem setting. This assumes to consider the equation (I4.38P in the extended space 

(4.42) (/, a) 6 X = #J x {a 6 M} and Q : X -> X, 

where proving the latter mapping for some compact subsets becomes a hard open problem. 
Note that here even the necessary convexity issue for applying Schauder's Theorem can 
be hard. We still do not know whether the representation such as (I4.42p may lead to any 
rigorous treatment of the nonlinear eigenvalue problem (14.381) for k > 1. 

5. General Cauchy problem: a homotopic approach 

5.1. Key concepts to justify: a first discussion. We now discuss some related prop- 
erties of more general solutions of the CP for the TFE (11.11) using a homotopic approach 
when the parameter n approaches zero. As shown in Section [31 we already know the 
similarity expression for the solutions of the "limiting" bi-harmonic equation (II .3p . This 
fundamental solution has also a self-similar structure thanks to the scaling invariance and 
the uniqueness of the fundamental solution of the equation (11.31) . denoted by (13.21) . 

The idea is to perform a homotopic approach from (11.11) to (II. 3p in order to reveal 
important (and still obscure in general) properties of the Cauchy problem. The reason 
is that the bi-harmonic equation (11.31) i.e., (II. ip when n — 0, with the same initial data, 
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admits the unique classic solution given by the convolution (13.11) . where b(x,t) is the 
fundamental solution (13.21) of the operator + A 2 , and the oscillatory rescaled kernel 
F(y) is the unique solution of the problem (13.31) . Hence, we expect that the knowledge 
of the solutions of ( 11. 3D can be extended to ( II. lj) at least for sufficiently small n > 0. In 
other words, we claim that the "fundamental" solutions for n = and small n > exhibit 
several similar properties, excluding, on the other hand, some others such as the compact 
support one for n > 0. In addition, the homotopic path n — > + can be used for a proper 
definition of the solutions of the Cauchy problem for the TFE-4 (II. ip . 

Thus, we assume that n > is sufficiently small. We define some "homotopic classes" 
of degenerate parabolic PDE's saying that the TFE (11.11) is homotopic to the linear PDE 
(II. 3p if there exists a family of uniformly parabolic equations (a homotopic deformation) 
with coefficient <fi £ {u) analytic in both variables u G M and e G (0, 1], 



We should point out that such a limit for nonnegative and not changing sign solutions, 
with various non-analytic (and non-smooth) regularizations has been widely used before 
in TFE-FBP theory as a key foundation; cf. ||, j|, and [7]. 
A possible homotopic path can be 



For any e G (0, 1], denote by u e (x, t) the unique solution of the CP for the regularized non- 
degenerate equation (15. ip with same data u . By classic parabolic theory, u £ is continuous 
(and analytic) in e G (0, 1] in any natural functional topology. The main problem is the 
behaviour as e — > 0, where the regularized PDE loses its uniform parabolicity. For second- 
order parabolic equations obeying the Maximum Principle, such regularization-continuity 
approaches are typical for constructing unique solutions with singularities (finite time 
blow-up, extinction, finite interfaces, etc.); see [18J as a source of key references and basic 
results. However, for higher-order degenerate parabolic flows admitting strongly oscilla- 
tory solutions of changing sign, such a homotopy-continuity approach generates a number 
of difficult problems. In fact, despite the fact that the passage to the limit as e — > looks 
like a reasonable way to define a proper solution of the TFE, we expect that there are 
always special classes of compactly supported initial data, for which such a limit is non- 
existent and, moreover, there are many partial limits, thus defining a variety of different 
solutions (meaning nonuniqueness) , as we show below. 

5.2. Preliminary estimates. To ascertain such a limit for (15. ip when e — > 0, we firstly 
obtain some estimations for its regularized solutions {u £ (x,t)}. Here, by Q we denote 
either WL N , or, equivalently, the bounded domain r D {t}, i.e., the section of the support. 




u\ n as £ — > uniformly on compact subsets. 



u t = -V • (0 £ (m)VAm) , 



e n + {l-e){e 2 + u 2 f^ eG (0,1]. 
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Proposition 5.1. Let u e (x,t) be the unique global solution of the CP for the regularized 
nondegenerate equation ( 15.11) with the initial data uq. Then, for any t G [0, T\, the 
following is satisfied: 

(i) u e (-,t) eH^Q); 

(ii) u £ (-,t) G L P {VL), with p = 1, 2, oo; and 

(iii) h £ G L 2 (tt x [0,T]), with h £ : = <j) £ (u £ )V Au £ . 

Proof. Firstly, multiplying (15.11) by Au e , integrating in Q x [0, t] for any t G [0, T], and 
applying the formula of integration by parts yield 

(5.3) | / \Vu £ (x,t)\ 2 + } j <P £ {u)\V Au £ \ 2 = \ j \Vu £ (x,0)\ 2 , 

n on n 

thanks to the boundary conditions (12.11) . Note that 

j[\Vu £ (x,t + h)\ 2 -\Vu £ (x,t)\ 2 ] 
n 

= — f[Au e (x, t + h) + Au £ (x, t)][u e (x, t + h) — u e (x, t)]. 

Then, dividing that equality by h, passing to the limit as h I 0, and integrating between 
and any t G [0, T], we find that 

/ / Au e u t = \ J \Vu e (x, t)\ 2 -\j \Vu e (x, 0)| 2 , 
on n n 

which provides us with the necessary expression to obtain (15.31) . Thus, from (15.31) . we have 
that (in fact, this is true from the beginning for classic C°°-smooth solutions of (15. ip . but 
we will need those manipulations in what follows) 

t 

(5.4) J\Wu £ (x,t)\ 2 < K and J J <f) £ (u)\V Au £ \ 2 < K, 
n on 

since both terms of the left-hand side in (15. 3p are always positive and the right-hand side 
is bounded by (12. 2p . for some positive constant K > that is independent of e. Then, 

Vu e (-,t) G L 2 (Q) for any t G [0,T]. 

Moreover, by Poincare's inequality, u £ (-,t) G L 2 (Q), and hence, 

(5.5) u e {-,t) G iZj}(fi) for any te[0,T). 
In fact, we may assume that 

(5.6) u e (-,t) e L°°(Q) for all t E (0,T). 

Indeed, for N = 1, this follows from (15.41) by Sobolev's embedding. For N > 2, this is 
a natural assumption inherited from the smooth analytic parabolic flow (15.11) . though its 
full proof sometimes can be a difficult issue; we refer to scaling and other techniques that 
may be convenient, [23j. From the conservation of mass assumption, we can also assure 
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that u £ (-,t) G L 1 (fi). Also, expression (15.31) combined with u e (-,t) G provides us 

with the estimate 

(5.7) /i e G L 2 (n x [0,T]) (/i E = £ ( M£ )VA Me ). 
Indeed, from (15.41) we find that 

/ f[e n + (1 - e){e 2 + w 2 )f ]|VAw e | 2 < so that 
o n 

(5.8) e n jj\VAu £ \ 2 < K and j j(e 2 + u 2 )%\VAu £ \ 2 < K, 

on on 

since £ G (0, 1) with a constant K > independent of e. Now, using Holder's inequality, 

}j\h £ \ 2 < 2e 2n j j \VAu £ \ 2 + 2jj(e 2 + u 2 )^(e 2 + u 2 )^\VAu £ \ 2 , 
on on on 

by dSSD and flU} (note also that u £ (-,t) G L x (fi)), we obtain (j5T7) . □ 
Furthermore, the following estimates are also ascertained: 

Lemma 5.1. Let u e (x, t) t/ie unique global solution of the CP for the regularized uniformly 
parabolic equation (15.11) with the initial data uq. Then, there exists some positive constant 
K > such that, for x\,x 2 G Q, independently of e and t, 

\u £ (xi,t) — u £ (x2,t)\ < K\xi — x 2 \2 for odd N, and 

\u e (xi,t) — u £ (x 2 , t)\ < K\xi — x 2 \ for even N. 

Proof. Thanks to Sobolev's inequality, we have that, for the odd dimension N > 1, 

||«e('»*)|| c? o,J (n j < C \\ U e(^ t )\\ H ^-iN {ny With ji = 0,j 3 = l,js = 2, • • • , 

for a positive constant C > 0. On the other hand, when the dimension N > 1 is even, 

IM'>*)||co,i(n) - C \\ U e('> ^11^-^(0)' With I" 2 = = 1 '^ 6 = 2 '"' ' 

for some positive constant C > 0. Hence, since Cg°(fi) is dense in W fe,2 (O) = # fe (ft) for 
k > 1, by the analytic smoothness of the solutions of the uniformly parabolic PDE (15.11) 
with analytic coefficients, for some positive constant K independently of e and t, 

\u £ (xi,t) — u £ (x 2 ,t)\ < K\x\ — x 2 \? for odd N, and 
\u £ (xi,t) — u £ (x 2 ,t)\ < K\x\ — x 2 \ for even N. □ 

Moreover, as was noted in [31 [5] and [7], u £ (x,t) is also Holder continuous in time with 
exponent ~, i.e., w £ (;t) G C°-8([0,T]) for the one dimensional case. However, we provide 
a new version for the ^-dimensional case. 
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Lemma 5.2. Let u e (x,t) be the unique global solution of the CP for the regularized equa- 
tion ( 15. II) with the initial data uq. Then, there exists a positive constant r > such 
that 

N 

\u £ {x, t 2 ) - u £ {x, tx) I < r \t 2 - t x I w> , 
if the dimension N is odd, and 

N 



\U £ (x,t 2 ) - U E (x,ti)\ < T\t 2 - ti 



3AT+2 



if the dimension N is even, for any t\,t 2 G [0,T], independently of e and x. 

Proof. First, consider a non-negative cut-off function (p G Cq°(M. n ) such that supp (<p) C Q 
and j Q <f — 1. Subsequently, multiplying 05. ip by a test function yj 7 G Qj^IR^), where 

with some xq G ^ and a constant 7 > to be properly chosen later on, integrating over 
Q x t 2 ] and applying the formula of integration by parts, we find that 



(5.9) 



*2 *2 

- / / y? 7 w e , t - / / Vy? 7 • (0 £ (m £ ) VAw e ) = 0, 
ti n ti n 

where - / / ^ 7 « £ ,t = / ^7(^(^2) - w e (*i))- 



On the other hand, by the present choice of y? 7 , we know that J Q <£> 7 = 1, since the 
Jacobian of is -. Then, we have that 

u e {x 0> t 3 )-u e (xo i ti) = J ip^(x)(u £ (x ,t 2 ) - u £ (x ,ti)) 

n 

= J v9 7 (x)(m £ (x , t 2 ) - u £ (x, t 2 ) + u £ (x, t 2 ) - u £ (x, t{) + u £ (x, ti) - u £ (x , ti)) 
n 

< f <p 7 (x)\u £ (x,t 2 ) - u £ (x ,t 2 )\ + I / <f 7 (x)(u £ (x,t 2 ) - u £ (x,ti))\ 
n n 

+ f (p 7 (x))\u £ (x ,ti) - u £ (x,tx)\. 
n 

Owing to Lemma I5TT1 and (15.91) on the time interval (ti, t 2 ), taking into account in the last 
inequality, we obtain 



*2 



\u £ (x ,t 2 ) - u £ (x ,ti)\ < 2K\x - x \2 + j j V<y5 7 • (4> £ (u)VAu £ 



\u £ (xo,t 2 ) - u £ (x ,ti)\ < 2K\x - x \ + 



;«)vam e ) 



with N odd, 



with N even. 



Moreover, by Holder's inequality and the choice of <£> 7 , 

ffV<p y -{<f> e (u)VAu e ) < (//|0 e (w)VAw e 
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i / *2 



ti n 



and, hence, thanks also to Proposition 15 .![ 



t 2 Jv+2 

// Vy? 7 • (0 £ (m)VAm £ ) < Tj-^*r\t 2 - h\2. 

h n 



Therefore, overall, interchanging x and Xq yields 



JV+2 



\u £ (x, t 2 ) — u £ (x, ti)\ < 2K\x — x \ 2 + T7 2N \t 2 — ti| 2 , with N odd, 

JV+2 1 

\u £ (x, t 2 ) — u £ (x, ti)\ < 2K\x — xo\ + r7 2N \t 2 — ti | 2 , with A" even. 
Thus, taking IK < r, |x — 3?o| < 7, and 7 < \t 2 — t\f, we obtain that 

|w e (a;, t 2 ) — w e (a;, ti)| < r\t 2 — ti|a +r|t 2 — *i| 2 ^ + ^, with iV odd, 

R (JV+2)/3 1 

|it E (x, ^2) — *i) I < r\t 2 — ti\ p + T\t 2 - ti\ 2N +2 , with even. 
Consequently, taking (3 = 2 (n+i) ' ^ ^ ^ s oc ^^' anc ^ ^ = 3^+2 wnen ^ i s even, completes 



the proof. □ 

5.3. Passing to the limit. To conclude this section, we show existence of weak solutions 
for the degenerate parabolic problem (11.11) passing to the limit as e goes to zero. However, 
we must admit from the beginning that, from the analysis performed below, it is not 
possible to assure which the limit will be (the solution of some CP or maybe the solution 
of some FBP). 

By Proposition EIU since, for bounded supports Q, the embedding Hq(Q) L 2 (Q) is 
compact, we can extract a convergent subsequence in L 2 (Q) as e [ for the solutions of 
(15. ip labelled again u £ (x,t) such that 

(5.10) lim || u e (-,t)- U(-,t)\\ L2{n) =0. 

Consequently, the convergence of the non-degenerate solutions of the problem (15. ip is 
strong in L 2 (Q). 

Moreover, thanks to the Holder continuity proved in Lemmas 15. II and 15.21 we have a 

- JV" — q -r JV — 

strong convergence as e | in C ' 2 ' 4 ( Ar + 1 ) (Q x [0, T]), when iV is odd, and in C ' ' 3JV + 2 (Q x 
[0,T]), when iV is even. This is possible after applying the Ascoli-Arzela Theorem, since 
{u e } is uniformly bounded and equicontinuous in Q x [0, T]. Of course, these estimates 
can imply other even stronger convergence results, which are not treated below in detail. 

Note that one difficulty we face is whether this limit depends on the taken subsequence 
or not. In other words, this analysis does not include any uniqueness result, which is 
expected to be a more difficult open problem for such nonlinear degenerate parabolic 
TFEs in non-fully divergence form and with non-monotone operators. However, the 
principal issue of the analytic regularization via (15.11) is that it is expected to lead to a 
smoother solution at the interface than those for the standard FBP. The difference is that 
the analytic regularized family {u £ }, in addition to (12. ip . is assumed to guarantee that, 
a.e. on the interface (assumed now sufficiently smooth), 

(5-11) = 0. 
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In fact, proper oscillatory solutions of the CP are assumed to exhibit even more regularity 
at smooth interfaces 



(5.12) = 0, where /=[£]- 1. 

Therefore, asn->0 + , the smoothness of such solutions at the interfaces increases without 
bounds. Obviously, this is not the case for the FBP (a "positive obstacle" one) with 
standard conditions (12. ip and a usual quadratic ( "parabolic" ) decay at the interfaces. 

Thus, as above and customary, multiplying (15. ip by a test function ip e C£°(Cl x (0, T)) 
and integrating by parts in Q x [0, T] gives 

T T 

-J J <p t u £ -JJVip- ((j) £ (u)VAu £ ) = 0. 
on on 

Next, operating with this equality, we find that 

T T T 

(5.13) / / <p t ue + e n J JV(p- VAu £ + (1 - e) J fV<p- ((e 2 + u 2 )^V Au £ ) = 0. 
o n on on 

Applying Holder's inequality, it is clear from (15.31) that there exists a subsequence labeled 
by {ek} such that the second term of (15. 13ft approximates zero as Ek | for a sufficiently 
small, n ~ 0, n > 0, 



eljjVp- VAu £k < e k (ef^ j j(VAu £k ) 2 ) ~ 2 ( / / | Vy,| 2 ) 1 < Ke k [ 0, 
on y on /v on 7 

as Ek I 0, for some positive constant K > 0. 

Moreover, on the subset V := {(x, t) 6 ft X [0, T] ; \u(x, t)\ > 5 > 0}, for any arbitrarily 
small 5 > 0, it is clear that the limiting solution as n — > is a weak solution of the 
bi-harmonic equation (11.31) . Indeed, by the regularity of the uniformly parabolic equation 
(15. ip and the uniformly Holder continuity of its solutions proved in Lemmas 15.11 and 15.21 
we obtain that u £tt , Vw e , Au £ , VAu £tX , and A 2 u £ converge uniformly on compact subsets 
of V . In general, it is not that difficult to see that, as e = — > (along the lines of 
classic results in [5] and related others), we obtain a weak solution of the TFE-4, i.e., 

(5.14) JfvtU + ffVtp- \U\ n VAU = 0, 

op op 

where U(x, t) is the limit obtained through (15.101) . We naturally assume that (p G C£°(V). 

However, in the "bad" subset {\u\ < 5}, for any sufficiently small 5 > 0, we must take 
e > sufficiently small and depending on 5. Indeed, we take e such that < e < 5. Thus, 
applying Holder's inequality to the third term in (I5.13P over the subspace where \u\ < 5, 
we have that 

j J y V -{{l-s)(e 2 + u 2 £ )^VAu £ 

{\u\<8} 

<(l I |V^| 2 Wf / {l-e) 2 {e 2 + u 2 £ T\VAu £ 

{|«|<5} 7 v {|«|<<5} 
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Then, since p G C^(Q x (0, oo)) and s G (0, 1), we get 

J / V^((l-£)( £ 2 + u 2 )SVAu £ ) <c(j J (l-e)(e 2 + ul) n \VAu 



{|u|<5} x {|u|<<5} 

for some positive constant C > 0. Making use of the fact that |it| < S and by (15. 4p . we 
find that 

j J V^((l-£)(e 2 + M£ 2 )iVA M£ ) < 

{\u\<8} 

C(j J (l-e)( £ 2 + 5 2 )t(e 2 + M 2 )f|VA M£ 

V {[«]<*} 

and, hence, using (15. 8ft . 

(5.15) J / V^-((l-e)(£ 2 + u 2 )fVA M£ ) <C 1( 5t -det, 

{|„[<*} 

for some constant C% > and taking e ~ 5 sufficiently small. 

Finally, the estimate (I5.15P shows the actual rate of the limit as n — > 0, together with 
e — > 0, in the analytic approximating flow (15.11) to get in this limit weak (and hence classic 
by standard parabolic theory) solutions of the bi-harmonic equation (11.31) . Namely, one 
has to have that 

(5.16) n = n(e) — > such that — » =^ ra(e) 3> ase^ + . 

However, this is not the end of the problem: indeed, under the condition (I5.16P on 
the parameters, we definitely arrive at the limit e, n(e) — > to the weak solution of the 
bi-harmonic equation written in the following "mild" form: 

(5-17) /„ tptU + J r /„ ■ VAU = 0. 

This is not a full definition of weak solutions, since it assumes just a single integration 
by parts, so allows us also positive solutions of the "obstacle" FBP for (jl.3p with the 
corresponding conditions (12 .ip (with n — 0), which can be constructed by "singular" 
regularization as in [5]. 

Thus, unfortunately, our analysis still does not recognize the desired difference between 
oscillatory solutions of the CP and others (possibly positive ones) of the standard FBP and 
others that can be posed for the TFE-4 (11.11) . Nevertheless, this first step in a homotopy 
analysis declares useful estimates and bounds on the parameters of regularization such 
as (15.161) . which are absolutely necessary for passing to the limit to get sign changing 
solutions of the linear bi-harmonic flow. 

Then, the homotopy concept as a connection to the linear PDE (15. ip can describe the 
origin (at n = 0) of the oscillatory solutions of TFEs and hence establish a transition to 
the maximal regularity of the solutions of (II. ip . Indeed, inevitably, bearing in mind the 
oscillatory character of the kernel -F(|y|) of the fundamental solution, the proper solutions 
of the CP are going to be oscillatory near finite interfaces at least for small n > 0. 
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5.4. Final remark: towards a full homotopy approach. The main obstacle faced 
in obtaining a full exhibition of such a homotopic approach, when n — > + , is the study 
of difficult pointwise limits near interfaces, where key singularities (a kind of Riemann's 
problems) occur. In particular, as a clue, let us mention a proper topology, in which 
we can define the homotopy. This is clear when we transform (11.11) (here, for simplicity, 
avoiding e-regularization as in (15. ip ) into a perturbation of the bi- harmonic equation (II. 3ft : 

(5.18) u t = -A 2 u +g n (u), where g n (u) := V • ((1 - \u\ n )VAu) , 

which we write down as an integral equation using the compact resolvent (B — A/) -1 and 
the semigroup from Section [31 We then deal with the integral equation 

(5.19) u{t) = b(t) * u + j'l b(t - s) * g n (u(s)) ds, 

where b(t) is the fundamental solution (13. 2p . As usual, the integral form (15.191) allows 
us to weaken the necessary treatment of the perturbation g n (u), which is assumed to be 
small as n — > 0. However, this does not rule out the principal difficulty concerning such 
an unusual and very sensitive perturbation g n {u). 

Thus, the main open problem of the homotopy issues is as follows: under which 
"topology-functional-geometric" setting for admitted solutions u(x,t), 



(5.20) J * b(t — s) * g n (u(s)) ds — > a.e. as n — > 0' 



In particular, it is not difficult to see that, for sufficiently smooth functions u with a finite 
number of transversal zero surfaces uniformly in small n > 0, we have that 

g n {u) — > as n — > + . 

at least a.e., and in other natural (weighted) topologies associated with the operator B 
and/or others. This can even be true uniformly on compact subsets, if the differential 
operators in u are bounded on such special functions, whose regularity and the "geo- 
metric transversal structure" near zero surfaces well-correspond to the desired maximal 
regularity/structure that are generic for the TFE solutions. But such a detailed a priori 
information on solutions seems to be excessive and not required. Taking the full homo- 
topic approach as e , n — ► 0, we desperately need to understand the structure of the zero 
surfaces, because of the oscillatory behaviour of these solutions of changing sign. Away 
from small neighbourhoods of such zeros (zero curves or surfaces), there is no any essen- 
tial problem. The structural properties of zeros of solutions of the TFE-4 (11.11) and the 
TFE-6 have been discussed in [T4"t [16] , and the results therein inspire us with a certain 
optimism concerning the correctness of the general homotopy approach to the CP for the 
TFEs, though difficulties are far away from being properly settled in a general setting. 
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